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In this contribution we review and discuss several aspects of Cosmic Voids. Voids are a 
major component of the large scale distribution of matter and galaxies in the Universe. 
Their instrumental importance for understanding the emergence of the Cosmic Web 
is clear. Their relatively simple shape and structure makes them into useful tools for 
extracting the value of a variety cosmic parameters, possibly including even that of 
the influence of dark energy. Perhaps most promising and challenging is the issue of the 
galaxies found within their realm. Not only does the pristine environment of voids provide 
a promising testing ground for assessing the role of environment on the formation and 
evolution of galaxies, the dearth of dwarf galaxies may even represent a serious challenge 
to the standard view of cosmic structure formation. 
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1. Cosmic Depressions 

Voids have been known as a feature of galaxy surveys since the first surveys were 
compiled' 10 ' - 3 Voids are enormous regions with sizes in the range of 20 — 50/i 
Mpc that are practically devoid of any galaxy, usually roundish in shape and oc- 
cupying the major sha re of space in the Universe. Forming an essential ingredient 
of the Cosmic Web^^, they are surrounded by elongated filaments, sheetlike walls 
and dense compact clusters. Following the discovery by Ref. l60l and IBTI of the most 
dramatic specimen, the Bootes void, a hint of their cen tral position within a weblike 
arrangement came with the first Cf A redshift slice ^1. This view has been dramati- 
cally en dorsed and expanded by the redshift maps of the 2dFGRS and SDSS surveys 
| 17 | 104 | Thgy na ve established voids as an integral component of the Cosmic Web. 
The 2dFGRS maps and SDSS maps (see e.g. fig. Q]) are telling illustrations of the 
ubiquity and prominence of voids in the cosmic galaxy distribution. 

In a void-based description of the evolution of the cosmic matter distribution, 
voids mark the transition scale at which density perturbations have decoupled from 
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Fig. 1. SDSS is the largest and most systematic sky survey in the history of astronomy. It 
is a combination of a sky survey in 5 optical bands of 25% of the celestial (northern) sphere. 
Each image is recorded on CCDs in these 5 bands. On the basis of the images/colours and their 
brightness a million galaxies are subsequently selected for spectroscopic follow-up. The total sky 
area covered by SDSS is 8452 square degrees. Objects will be recorded to m,[ im = 23.1. In total 
the resulting atlas will contain 10 s stars, 10 s galaxies and 10 5 quasars. Spectra are taken of 
around 10 6 galaxies, 10 5 quasars and 10 5 unusual stars (in our Galaxy). Of the 5 public data 
releases 4 have been accomplished, ie. 6670 square degrees of images is publicly available, along 
with 806,400 spectra. In total, the sky survey is now completely done (107%), the spectroscopic 
survey for 68%. This image is taken from a movie made by Subbarao, Surendran & Landsberg 
(see website: http://astro.uchicago.edu/cosmus/projects/sloangalaxies/). It depicts the resulting 
redshift distribution after the 3rd public data release. It concerns 5282 square degrees and contained 
528,640 spectra, of which 374,767 galaxies. 



I 
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the Hubble flow and contracted into recognizable structural features. On the basis 
of theoretical models of void formation one might infer that voids may act as the key 
organizin g element fo r arranging matter concentrations into an all-pervasive cosmic 
network | 55 | 86 | 114 | 97 | ^ g vo j(j s expand, matter is squeezed in between them, and 
sheets and filaments form the void boundaries. This view is supported by numerical 
studies and computer simulations of the gravitational evo lution of voids in more 
complex and realistic configurations | 68 | 86 | 24 | 115 | 33 | 16 | 71 | mar ]j ec j example of 
the evolution of a typical large and deep void in a ACDM scenarios is given by the 
time sequence of six frames in fig. [5J 

Soon after their discovery, various studies poi nted o ut their essential role in 
the organization of the cosmic matter distribution I 55 | 86 [ xheir effective repulsive 
influence over their surroundings has been recognized in various galaxy surveys in 
the Local Universe. 



1.1. Voids: cosmological importance 

There are a variety of reasons why the study of voids is interesting for our under- 
standing of the cosmos. 

• They are a prominent aspect of the Megaparsec Universe, instrumental in the 
spatial organization of the Cosmic Web. 

• Voids contain a considerable amount of information on the underlying cosmo- 
logical scenario and on global cosmological parameters. Notable cosmological 
impri nts are fo und in the outflow velocities and accompanying redshift distor- 
tions HUSHED], Also their intrinsic structure, shape and mutual a lignment arc 
sensitive to the cosmology, including that of dark energy [ Z4 | 66 | 81 | -j^g cogmo _ 
logical ramifications of the reality of a supersized voids akin to the identified 
ones by ref. [89] and ref. [35] would obviously be far-reaching. 

• The pristine low-density environment of voids represents an ideal and pure set- 
ting for the study of galaxy formation and the influence of cosmic environment 
on the formation of galaxies. Voids are in particular interesing following the ob- 
servation by Peebles that the dearth of low luminosity objects in voids is hard 
to understand within the ACDM cosmology ^1 



1.2. Void Sizes 

For the most systematic and complete impression of the cosmic void population 
the Local Universe provides the most accessible region. Recently, the deep view 
of the 2dFGRS and SDSS probes (see fig. [IJ has been supplemented with high- 
resolution studies of voids in the nearby Universe. Based upon the 6dF survey ^Sl, 
Fairall (person, commun.) identified nearly all void s w ithin the surveyed region out 
to 35, 000 km s _1 . It is the 2MASS redshift survey^- the densest all-sky redshift 
survey available - which has provided a uniquely detailed and complete census of 
large scale structures in our Local Universe 
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Voids in the ga laxy distribution account for about 95% of the total volume 
| 58 | 27 | 28 | 49 | 83 | 88 | 8 0] The typical sizes of voids in the galaxy distribution depend on 
the galaxy population used to define the voids. Voids defined by galaxies brighter 
than a typical L* galaxy tend to have diameters of order 10 — 20/i~ 1 Mpc, but voids 
associated with rare luminous galaxies can be con sidera bly larger; diameters in the 
range of 20/i _1 — 50/i _1 Mpc are not uncommon | 49 | 8 3] pj rm U pp 6 r limits on the 
maximum void size have not yet been set. Recently there have been claims of the 
existence of a supersized void, in the counts of the NWS catalogue of radio sourcs, 
and of its possible imprint on the CMB via the ISW effect in the form of a 5° CMB 
'Cold Spot' SS. A systematic search for such supervoids, and superclusters, in the 
Luminous Red Galaxy (LRGs) sample from the SDSS survey has indeed provided 
one of the strongest claims for a significant (4cr) detection of the ISW effect in the 
WMAP maps I 35 | 36 |_ jf ^iis will be confirmed it will pose an interesting challenge 
to concordance cosmological scenarios ^1. 

At the low end side of the void size distribution, a few studies claim to have found 
what may be the smallest _genuine voids in existence. On the basis of the Catalog 
of Neighbouring Galaxies ref. 11061 identified a total of some 30 minivoids, each 
completely devoid of galaxies, with sizes of 0.7 — 3.5/i _1 Mpc. 



2. Formation and Evolution of Voids 

At any cosmic epoch the voids that dominate the spatial matter distribution are a 
manifestation of the cosmic structure formation process reaching a non-linear stage 
of evolution. 

Voids emerge out of the density troughs in the primordial Gaussian field of den- 
sity fluctuations. Early the oretical m odels of void formation concentrated on the 
evolution of isolated voids | 47 | 55 | 5 | 8 ] a resm |; G f their underdensity voids rep- 
resent a region of weaker gravity, resulting in an effective repulsive peculiar grav- 
itational influence. Initially underdense regions therefore expand faster than the 
Hubble flow, and thus expand with respect to the background Universe. Fig. [^illus- 
trates the evolution of a spherical isolated void. As matter streams out of the void, 
the density within the void decreases, with isolated voids asymptotically evolving 
towards an underdensity 8 = — 1, pure emptiness. The same expanding and evac- 
uating behaviour of void regions apply in the far more complex circumstances of 
the real cosmic matter distribution. The illustration of a void in a ACDM Universe 
is illustrated in fig. [2] by a sequence of 6 timesteps. Because the density within 
underdense regions gradually increases outward, we see a decrease of the corre- 
sponding peculiar (outward) gravitational acceleration: void matter in the centre 
moves outward faster than void matter towards the boundary. This leads to matter 
accumulating in ridges which surround the void, while the interior evolves into a 
uniform low-density region resembling a low-density homogeneous FRW Universe. 
The steepness of the resulting density profile depends on the protovoid depression 
^1 In nearly all conceivable situations the void therefore appears to assume a bucket 
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shape, with a uniform interior density depression and a steep outer boundary (fig. [2J 
right frame). 

A characteristic evolutionary timescale for voids is that of shellcrossing. This 
happens when interior shells of matter take over initially exterior shells. Ref. [5| 
demonstrated that once voids have passed the stage of shellcrossing they enter 
a phase of self-similar expansion (fig. [2]). Subsequently, their expansion will slow 
down with respect to the earlier linear expansion. This impelled ref. [8] to identify 
voids in the present-day galaxy distribution with voids that have just reached the 
stage of shell-crossing. It happens when a primordial density depression attains 
a linearly extrapolated underdensity S v = f v = —2.81 (for an EdS universe, see 
sec !4.1|) . A perfectly spherical "bucket" void will have expanded by a factor of 1.72 
at shellcrossing, and therefore have evolved into an underdensity of ~ 20% of the 
global cosmological density, ie. S v>n i = —0.8. In other words, the voids that we 
see nowadays in the galaxy distribution do probably correspond to regions whose 
density is ~ 20% of the mean cosmic density (note that it may be different for 
underdensity in the galaxy distribution). 

Note that while by definition voids correspond to density perturbations of at 
most unity, \S V \ < 1, mature voids in the nonlinear matter distribution do repre- 
sent highly nonlinear features. This may be best understood within the context of 
Lagrangian perturbation theory ^2. Overdense fluctuations may be described as a 
converging series of higher order perturbations. The equivalent perturbation series 




Fig. 2. Spherical model for the evolution of voids. Left: a pure (uncompensated) tophat void 
evolving up to the epoch of shell-crossing. Initial (linearly extrapolated) density deficit was 
^Un.o = —10.0, initial (comoving) radius R, o = 5.0h~ 1 Mpc. Right: a void with an angular 
averaged SCDM profile. Initial density deficit and characteristic radius are same as for the tophat 
void (left). The tendency of this void to evolve into a tophat configuration by the time of shell 
crossing is clear. Shell-crossing, and the formation of a ridge, happens only if the initial profile is 
sufficiently steep. 
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Fig. 3. Simulation of evolving void (LCDM scenario). A void in a n = power-law power spectrum 
model. The slice is 50h~ x Mpc wide and 10h _1 Mpc thick. Shown are the partciles and smoothed 
density field (smoothed on a scale of 4/i _1 Mpc) at six different timesteps: a=0.05, 0.15, 0.35, 0.55, 
0.75 and 1.0. Image courtesy of Erwin Platen 
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is less well behaved for voids: successive higher order terms of both density deficit 
and corresponding velocity divergence alternate between negative and positive. 

3. Void Dynamics 

Figure Q] shows a typical void-like reg ion in a ACDM Universe. It concerns a 256 3 
particles GIF iV-body simulation I^HI encompassing a ACDM (f2 TO = 0.3, f^A = 
0.7, Hq — 70 km/s/Mpc) density field within a (periodic) cubic box with length 
141/i _1 Mpc and produced by means of an adaptive P 3 M A^-body code. 

The top left frame shows the particle distribution in and around the void within 
this 42.5/i _1 Mpc wide and l/i _1 Mpc thick slice through the simulation box. In the 
same figure we include panels of the density and velocity field in the void, determined 
by means of a DTFE reconstruction E 5 l 93 | 116 | -j-he region appears as a slowly 
varying region of low density (top righthand frame). Notice the clear distinction 
between the empty(dark) interior regions of the void and its edges. In the interior 
of the void several smaller subvoids can be distinguished, with boundaries consisting 
of low density filamentary or planar structures. 

The general characteristics of the expanding void are most evident when fol- 
lowing the density and velocity profile along a one-dimensional section through the 
void. The bottom- left frame of fig. 2] shows these profiles for the linear section along 
the solid line indicated in the other three frames. The density profile does confirm 
to the general trend of low-density regions to develop a near uniform interior den- 
sity surrounded by sharply defined boundaries. Nonetheless, we see that the void is 
interspersed by a rather pronounced density enhancement near its centre. 

The flow in and around the void is dominated by the outflow of matter from the 
void, culminating into the void's own expansion near the outer edge. The compari- 
son with the top two frames demonstrates the strong relation with features in the 
particle distribution and the density field. Not only it is slightly elongated along 
the direction of the void's shape, but it is also sensitive to some prominent internal 
features of the void. Towards the "SE" direction the flow appears to slow down near 
a ridge, near the centre the DTFE reconstruction identifies two expansion centres. 

The void velocity field profile is intimately coupled to that of its density field. The 
linear velocity increase is a manifestation of its general expansion. The near constant 
velocity divergence within the void conforms to the super-Hubble flow expected for 
the near uniform interior density distribution. Because voids are emptier than the 
rest of the universe they will expand faster than the rest of the universe with a net 
velocity divergence equal to 







H 



3(a- 1), 



a — Hvoid/ H , 



(1) 

(2) 



where a is defined to be the ratio of the super-Hubble expansion rate of the void 
and the Hubble expansion of the universe. 
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Fig. 4. The density and velocity field around a void in the GIF LCDM simulation. The top 
righthand panel shows the N-body simulation particle distribution within a slice through the 
simulation box, centered on the void. The top righthand panel shows the grayscale map of the 
DTFE density field reconstruction in and around the void, the corresponding velocity vector plot 
is shown in the bottom lefthand panel. Notice the detailed view of the velocity field: within the 
almost spherical global outflow of the void features can be recognized that can be identified with 
the diluted substructure within the void. Along the solid line in these panels we determined the 
linear DTFE density and velocity profile (bottom righthand frame). We can recognize the global 
"bucket" shaped density profile of the void, be it marked by substantial density enhancements. 
The velocity field reflects the density profile in detail, dominated by a global super-Hubble outflow. 
From Schaap 2007. 



3.1. Dynamical Influence 

Various studies have found strong indications for the imprint of voids in the peculiar 
velocity flows of galaxies in the Local Universe. Ref. 13. made the first claim of seeing 
pushing influence of voids when assessing the stronger velocity flows of galaxies 
along a filament in the first CfA slice. Stronger evidence came from th e ext ensive 
and systematic POTENT analysis of Mark III peculiar galaxy velocities in the 



I 
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Fig. 5. Gravitational impact of the Sculptor Void. The righthand frame shows the inferred velocity 
field in and around the Sculptor void near the Local Supercluster. The colour map represents 
the density values, with dark blue at <5 ~ —0.75 and cyan near <5 ~ 0.0. The vectors show the 
implied velocity flow around the void, with a distinct nearly spherically symmetric outflow. It 
is a zoom-in onto the indicated region in the density and velocity map in the Local Universe 
(lefthand) determined on the basis of the PSCz galaxy redshift survey. From: Romano-Diaz & van 
de Weygaert 2007. 



Local Universe EHEJ POTENT found that for a fully selfconsistent reconstruction 
of the dynamics in the Local Universe, it was inescapable to include the dynamical 
influence of voids 

With the arrival of new and considerably improved data samples the dynamical 
influence of voids in the Local Universe has been investigated and understood in 
greater detail. The reconstruction of the density and velocity field in our local cosmos 
on the basis of the 2MASS redshift survey has indeed resulted in a very interesting 
and complete view of the dynamics on Megaparsec scales. This conclusion agree 
with that reached on the basis of an analysis of the peculiar velocity of the Local 
Group by ref. 11091 Their claim is that the Local Void is responsible for a considerable 
repulsive influence, accounting for ~ 259 km s _1 of the ~ 631 km s _1 Local Group 
motion with respect to the CMB. 
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Fig. 6. Voids and Void Dynamics in a ACDM cosmology. In a 100/i - 1 Mpc box, the lh~ 1 Mpc 
Gaussian smoothed density field contour levels are superimposed on the lh — 1 Mpc smoothed map of 
the tidal force field strength (note the outstanding cluster locations), along with ellipses identifying 
the voids by size, shape and orientation (on a scale of 6h~ 1 Mpc). 



3.2. Void Shapes and Tides 

Ref. [5S] pointed out that any (isolated) asphcrical underdensity will become more 
spherical as it expands. The effective gravitational acceleration is stronger along the 
short axis than along the longer axes. For overdensities this results in a stronger 
inward acceleration and infall, producing increasingly flattened and elongated fea- 
tures. By contrast, for voids this translates into a larger outward acceleration along 
the shortest axis so that asphericities will tend to diminis h. Fo r the interior of voids 
this tendency has been confirmed by N-body simulations Hl^l 

In reality, voids will never reach sphericity. Even though voids tend to be less 
flattened or elongated than the halos in the dark matter distribution, they are 
nevertheless quite nonspherical: ref. [81] find that they are slightly prolate with axis 
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ratios of the order of c : b : a w 0.5 : 0.7 : 1. This agrees with the statement by ref.[95l 
and ref. [74] that in realistic cosmological circumstances voids will be nonspherical. 
This is also quite apparent in images of, for example, the Millennium simulation ESI. 

The flattening is a result of large scale dynamical and environmental factors, 
amongst which we can identify at least two important factors Even while their 
internal dynamics pushes them to a more spherical shape they will never be able 
to reach perfect sphericity before encountering surrounding structures such as over- 
dense filaments or planar walls. Even more important may be the fact that, for 
voids, external tidal influences remain important: voids will always be rather mod- 
erate densities perturbations since they can never be more underdense than 5 = — 1. 
External tidal forces are responsible for a significant anisotropic effect in the devel- 
opment of the voids. In extreme cases they may even cause complete collapse of the 
void. 




0.5 1 1.5 2 2.5 0.5 1 1.5 2 2.5 
a(t) a(t) 

Fig. 7. Anisotropic Void Expansion and Collapse. Evolution of the same initially spherical void. 
Dashed line: without external influences Solid line: under the influence of an axisymmetric and 
linearly evolving external tidal force field (-E mmj o = ( — E, —E, 2E)). Topleft: (comoving) void axis, 
c' m = c m /a(t). Topright: axis ratios C2/C1 and C3/C1. Bottom Left: Internal density p e of the void, 
in units of the cosmic density p u (t). Bottom Right: The velocity v m along the axes of the voids, 
in units of the Hubble velocity vh- Note the collapse along axis 1 and 2. 
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3.2.1. Ellipsoidal Voids 



The major influence of external tidal forces on the evolution of voids can be un- 
derstood by assessing the evolution of homoge neous underdens e ellipsoids ^1 To 
some extent, the Homogeneous Ellipsoidal represents a better 

description for the evolution of voids than that of collapsing overdense peaks. Over- 
dense regions contract into more compact and hence steeper density peaks, so that 
the area in which the ellipsoidal model represents a reasonable approximation will 
continuously shrink. On the other hand, we have seen the natural tendency of voids 
is to develop an increasingly uniform interior density field (see fig. [2]): while they 
expand their interior gets drained of matter and develops a flat "bucket-shaped" 
density profile. 

In many respects the homogeneous model is a better approximation for under- 
dense regions than it is for overdense ones. Evidently, the approximation is restricted 
to the interior and fails at the void's outer fringes because of its neglect of the dom- 
ineering role of surrounding material, such as the sweeping up of matter and the 
encounter with neighbouring features. 

The homogeneous ellipsoidal model assumes an object to be a region with a 
triaxially symmetric ellipsoidal geometry and a homogeneous interior density, em- 
bedded within a uniform background density p u . Consider the simple situation of 
the external tidal shear directed along the principal axes of the ellipsoid. The gravi- 
tational acceleration along the principal axes of an ellipsoid with over/underdensity 
5 can be evaluated from the expression for the corresponding scale factors 1Z U 

d 2 K u 



dr 



1 + 6 1 , 2 



Km ~ T m K m + AR m , (3) 



3 2 w 3 y 

where we have also taken into account the influence of the cosmological constant A. 
The factors a m (t) are the ellipsoidal coefficients specified by the integral equation, 

7^3 VR- ( 4 ) 

° a) n AKw+^r 

AJ -n— 1 

The influence of the external (large-scale) tidal shear tensor T^*' enters via the 
eigenvalue r m . From eqn. [3) it is straightforward to appreciate that as 6 grows 
strongly nonlinear, the relative influence of the large-scale (near-)linear tidal field 
will decline. However, the density deficit S of voids will never exceed unity, \S\ < 1, 
so that the importance of the factor r m remains relatively large. Its impact can be 
so strong that it not only effects an anisotropic expansion of the void, but even may 
manage to make an initially spherically expanding void to collapse (see fig. E]). 

3.2.2. Void Alignments 

Large scale tidal influences not only manifest themselves in the shaping of individ- 
ual voids. They are also responsible for a distinct alignment of substructures along 
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a preferred direction, while they are also instrumental in their mutual arrange- 
ment and organization. Locally, the orientation of a void turns out to be strongly 
aligned with the tidal force field generated by structures on scales up to at least 
20 — 30/i _1 Mpc. This goes along with a similar mutual alignment amongst voids 
themselves. They have strongly correlated orientations over distances > 30/i -1 Mpc 
a scale considerably exceeding the typical void size. It forms a strong confirma- 
tion of the large scale tidal force field as the dominant agent for the evolution and 
spatial organization of the Megaparsec Universe, as emphasized by the Cosmic Web 
theory L^M 



4. Void Sociology 

Computer simulations of the gravitational evolution of voids in realistic cosmological 
environments do show a considerably more complex situation tha n that described 
by idealized spherical or ellipsoidal models |68 | 86|24 | 115 | 33 | 16 | 71 | j n receir t; y ears 
the huge increase in computational resources has enabled N-body simulations to 
resolve in detail the intricate substructure of v oids within the con text of hierar- 
chical cosmological structure formation scenarios | 70 | 34 | 46 | 2 | 33 | 16 | 7 T] xhey confirm 
the theoretical expectation of voids having a rich substructure as a result of their 
hierarchical buildup (see e.g. fig. H]). 

Ref. 97. treated the emergence and evolution of voids within the context of hi- 
erarchical gravitational scenarios. It leads to a considerably modified view of the 
evolution of voids. The role of substructure within their interior and the interaction 
with their surroundings turn out to be essential aspects of the hierarchical evolution 
of the void population in the Universe. An important guideline are the heuristic void 
model simulations by ref. [24] and the theoretical void study by ref. [91] within the 
context of a Lagrangian adhesion model approach by ref. [9TJ As argued by ref. [97l 
the hierarchical development of voids, akin to the evo lution of overdense halos, may 
be described by an excursion set formulation I 85 | 9 J96| 

In some sense voids have a considerably more complex evolutionary path than 
overdense halos. Two processes dictate the evolution of voids: their merging into 
ever larger voids as well as the collapse and disappearance of small ones embedded 
in overdense regions (see fig. [8|) . 



4.0.3. Void Merging 

First, consider a small region which was less dense than the critical void density 
value. It may be that this region is embedded in a significantly larger underdense 
region which is also less dense than the critical density. Many small primordial den- 
sity troughs may exist within the larger void region. Once small density depressions 
located within a larger embedding underdensity have emerged as true voids at some 
earlier epoch, their expansion tends to slow down. 

When the adjacent subvoids meet up, the matter in between is squeezed in 
thin walls and filaments. The peculiar velocities perpendicular to the void walls 
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10 12 14 16 10 13 14 16 10 12 14 16 
h- 1 Mpc h- 1 Mpc h" 1 Mpc 

Fig. 8. The two modes of void evolution: void merging (top row) and void collapse (bottom 
row). Top: three timesteps of evolving void structure in a 128 3 particle N-body simulation of 
structure formation in an SCDM model (a cxp = 0.1,0.3,0.5). The sequence shows the gradual 
development of a large void of diameter fs 25h — 1 Mpc as the complex pattern of smaller voids and 
structures which had emerged within it at an earlier time, merge with one another. It illustrates 
the void-in-void process of the evolving void hierarchy. Bottom: a choice of three collapsing voids 
in a constrained N-body simulation, each embedded within an environment of different tidal shear 
strength. The arrows indicate the velocity vectors, showing the infall of outer regions onto the void 
region. As a result the voids will be crushed as the surrounding matter rains down on them. 



are mostly suppressed so that the flow of matter is mostly confined to tangential 
motions. The subsequent merging of the subvoids is marked by the gradual fading 
of these structures while matter evacuates along the walls and filaments towards 
the enclosing boundary of the emerging void (also see top row fig. [8]). The 
timescale on which the internal substructure of the encompassing void is erased 
is approximately the same as that on which it reaches maturity. 

The final result is the merging and absorption of the subvoids in the larger void 
emerging from the embedding underdensity. Hence, as far as the void population is 
concerned only the large void counts, while the smaller subvoids should be discarded 
as such. Only a faint and gradually fading imprint of their original outline remains 
as a reminder of the initial internal substructure. 



I 
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4.0.4. Void Collapse 

A second void process is responsible for the radical dissimilarity between void and 
halo populations. If a small scale minimum is embedded in a sufficiently high large 
scale density maximum, then the collapse of the larger surrounding region will 
eventually squeeze the underdense region it surrounds: the small-scale void will 
vanish when the region around it has collapsed completely. Alternatively, though 
usually coupled, they may collapse as a result of the tidal force field in which they 
find themselves. If the void within the contracting overdensity has been squeezed to 
vanishingly small size it should no longer be counted as a void (see fig. El bottom 
row). 

The collapse of small voids is an important aspect of the symmetry breaking 
between underdensities and overdensities. In the primorial Universe, Gaussian pri- 
mordial conditions involve a perfect symmetry between under- and overdense. Any 
inspection of a galaxy redshift map or an N-body simulation shows that there is 
a marked difference between matter clumps and voids. While the number density 
of halos is dominated by small objects, void collapse is responsible for the lack of 
small voids. 



4.1. Void Excursions 

Marked by the two processes of merging and collapse, the complex hierarchical 
buildu p of the void population may be modelled by a two-barrier excursion set for- 
malism The barriers refer to the critical (linear) density thresholds involved with 
the merging and collapse of voids. Whenever the linearly extrapolated (primordial) 
8l{t, t\R) on a scale R, 

S L (r,t\R) = jM-6 L (r,ti\R), (5) 

exceeds a density threshold f c it will collapse. For an Einstein-de Sitter fi TO = 1 
Universe the critical value has the well-known value f c ~ 1.686 SOI. A void will 
form when an underdensity reaches the critical density threshold of shell- crossing, 
corresponding to a value of f v ~ —2.81 for spherical voids in an Einstein-de Sitter 
Universe. The linear density gro wth factor D (t), normalized to unity at the present 
epoch, follows from the integral [ 44 | 77 | 42 |64]^ 

D(t) = D(t, n m . , n At0 ) = H(a) P . (6) 

£ Jo a tl (a ) 

Emerging from a primordial Gaussian random field, many small voids may coexist 
within one larger void. Small voids from an early epoch will merge with one another 
to form a larger void at a later epoch. The excursion set formalism takes account 
of this void-in-void configuration by discarding these small voids from the void 
count. To account for the impact of voids disappearing when embedded in collapsing 
regions, the two-barrier formalism also deals with the void-in-cloud problem. 
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By contrast, the evolution of overdensities is governed only by the cloud-in- 
cloud process; the cloud-in-void process is much less important, because clouds 
which condense in a large scale void are not torn apart as their parent void expands 
around them. This asymmetry between how the surrounding environment affects 
halo and void formation is incorporated into the excursion set approach by using 
one barrier to model halo formation and a second barrier to model void formation. 
Only the first barrier matters for halo formation, but both barriers play a role in 
determining the expected abundance of voids. 

This image of void formation in the dark matter distribution has been elaborated 
by ref. [30] to describe the implications for voids in the galaxy distribution. 

4.2. Voidpatch Theory 

The extended Press-Schechter formulation of hierarchical void evolution is a local 
description and fails to take into account the proper spatial and geometric setting 
and configuration underlying the process. 

The peakpatch formalism introduced by ref. [10] represents a more proper non- 
local description of the hierarchical buildup of structures. Framed in an elaborate 
analytical framework it involves a natural definition of a peaks physical extension, 
and takes into the account of the anisotropic nature of gravitational collapse and of 
the relative location and size of peaks in and around evolving density peaks. 

Along similar lines ref. [82] have phrased a voidpatch formalism to follow the 
hierarchical buildup of voids. It will allow a much better understanding of the con- 
figuration, shape and substructure of individual voids emerging in a cosmic density 
field. It will also help towards the identification of collapsing subvoids and help 
towards evaluating the impact of the various dynamical influences. 

4.3. Void Population Statistics. 

The analytical evaluation of the two-barrier random walk problem in the extended 
Press-Schechter approach leads directly to a prediction of the distribution function 
n v (M) for voids on a mass scale M (or corresponding void size . The resulting 
void spectrurrQ. is peaked, with a sharp cutoff at both small and large values of the 
peak mass M Vi * (fig. [9] bottom lefthand frame), 

n v {M)AM » 

(7) 

^^expp^ exp(-^4-2^1, 

a The conversion of the void mass scale to equivalent void radius R is done by assuming the simplest 
approximation, that of the spherical tophat model. According to this model a void has expanded 
by a factor of 1.7 by the time it has mature, so that V v = (M/p u ) * 1.7 3 . 

b Note that for near-empty voids the mass scale M is nearly equal to the corresponding void mass 
deficit. 
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Fig. 9. Distribution of void radii predicted by the two-barrier extended PS formalism, in an 
Einstein de-Sitter model with P(k) oc fc — 15 , normalized to = 0.9 at z = 0. Top left panel shows 
the mass fraction in voids of radius r. Bottom left panel shows the number density of voids of 
radius r. Note that the void-size distribution is well peaked about a characteristic size provided 
one accounts for the void-in-cloud process. Top right panel shows the cumulative distribution of 
the void volume fraction. Dashed and solid curves in the top panels and bottom left panel show the 
two natural choices for the importance of the void-in-cloud process discussed in the text: <5 C = 1.06 
and 1.686, with 8 V = —2.81. Dotted curve shows the result of ignoring the void-in-cloud process 
entirely. Clearly, the number of small voids decreases as the ratio of 8 C /\8 V \ decreases. Bottom right 
panel shows the evolution of the cumulative void volume fraction distribution. The three curves in 
this panel are for <5 C = 1.686(1 + z), where z = (solid), 0.5 (dotted) and 1 (dashed). From Sheth 
&; van de Weygaert 2004. 



where \J v v (M) is the fractional relative underdensity, 



v v (M) 



l/v 



a(M) ' 



(8) 



in which with the dependence on the mass scale M entering via the rms density 
fluctuation on that scale, a(M). The quantity T> is the 11 void- and- cloud parametei } \ 
V = \f v \/(f c + |/v I)- It parameterizes the impact of halo evolution on the evolving 
population of voids: the likelihood of smaller voids being crushed through the void- 
in-cloud process decreases as the relative value of the collapse barrier / c with respect 
to the void barrier / v becomes larger. 



I 



I 
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From fig. [9] we also see that the population of large voids is insensitive to the 
void-in-cloud process. The large mass cutoff of the void spectrum is similar to the 
ones for clusters and reflects the Gaussian nature of the fluctuation field from which 
the objects have condensed. The characteristic void size increases with time: the 
gradual merging of voids into ever larger ones is embodied in a self-similar shift of 
the peak of the void spectrum. 

When evaluating the corresponding fraction of contained in voids on mass scale 
M, f(M) = Mn v (M)/p u (top lefthand frame, hg. [9|), we see that this also peaks 
at the characteristic void scale. It implies a mass fraction in voids of approximately 
thirty percent of the mass in the Universe, with most of the void mass to be found 
in voids of this characteristic mass. 

While the two-barrier excursion set formalism offers an attractive theoretical 
explanation for the distinct asymmetry between clumps and voids and for the peaked 
void size distribution, realistic cosmological simulations are needed to identify where 
the disappearing small-scale voids are to be found in a genuine evolving cosmic 
matter distribution. Using the GIF N-body simulations of various CDM scenarios, 
ref. [79] has managed to trace various specimen of this unfortunate void population. 
Subvoids in the interior of a larger void tend to merge with surrounding peers while 
the ones near the boundary get squeezed out of existence. These do not collapse 
isotropically, but tend to get sheared by their surroundings. 

An important aspect of the implied void population is that it is approximately 
space-filling. It underlines the adagio that the large scale distribution of matter may 
be compared to a soapsud of expanding bubbles. This follows from evaluation of the 
cumulative integral 



where the factor 1.7 is an estimate of the excess expansion of the void based upon 
the spherical model for void evolution (see footnote). The top righthand panel of 
fig. O shows the resulting (current) cumulative void volume distribution: for a finite 
value of void radius R the whole of space indeed appears to be occupied by voids, 
while the bottom righthand frame shows the gradual shift of the cumulative volume 
distribution towards larger voids. In other words, the correct image appears to be 
that of a gradually unfolding bubbly universe in which the average size of the voids 
grows as small voids merge into ever larger ones. 

4.4. Void Substructure 

An important issue within the hierarchically proceeding evolution of voids and the 
Cosmic Web is the fate of its substructure. In voids the diluted and diminished 
infrastructure remains visible, at ever decreasing density contrast, as cinders of 
the earlier phases of the void hierarchy in which the substructure stood out more 
prominent 
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N-body simulations show that voids do retain a rich infrastructure. Examples 
such as the top row of fig. [8] an d images of the Millennium simulation ESI show that 
while void substructure does fade, it does not disappear. We may find structures 
ranging from filamentary and sheetlike structures to a population of low mass dark 
matter halos and galaxies (see e.g. ref. II 151 &; [34]) . Although challenging, these may 
also be seen in the observational reality. The galaxies that populate the voids do 
currently attract quite some attention (see next section). Also, the SDSS galaxy 
survey has uncovered a substantial level of substructure within the Bootes void, 
confirming tantative indications for a filamentary feature by ref. 11021 

The gradual fading of the initially rich substructure as a void expands, and its 
matter content evacuates its interior, is clearly visible in the development of the 
void region in fig. [2] This may be seen as the slowing of structure growth in the 
low-density environment of voids. In their analytical treatment, ref. [33] described 
how the evolving substructure in a void could be translated into structure evolving 
in a lower Universe, including the implicit change in spectral characteristics of 
the density field (see ref. [69] for a reappraisal and criticism of some of the technical 
issues) . 

The resulting observation is one in which the density distribution in a void 
region undergoes a manifest and continuous transformation towards an ever larger 
dominating scale, while it mostly retains the topological characteristics of the initial 
density field. In a sense, it reflects of the buildup of the Cosmic Web itself, its basic 
pattern already recognizable in the primordial matter distribution HH. 



5. Void Galaxies 

A major point of interest is that of the galaxies populating the voids, the void galax- 
ies. The pristine environment of voids represents an ideal and pure setting for the 
study of galaxy formation. Largely unaffected by the complexities and processes 
modifying galaxies in high-density environments, the isolated void regions must 
hold important clues to the formation and evolution of galaxies. This has made 
the study of the relation between void galaxies and their surroundings an impor- 

tant aspect of the recent interest in environmenta l influences on galaxy formation 
I102I63I84I56I38I39I49I50I87I88I105I75I76I14I113I100I 

Amongst the issues relevant for our understanding of galaxy and structure for- 
mation, void galaxies have posed several interesting riddles and questions. Arguably 
the most prominent issue is that of the near absence of low- luminosity galaxies in 
voids, while ACDM cosmology expects voids to be teeming with dwarfs and low sur- 
face brightness galaxies. A manifest environmental influence which still needs to be 
understood is the finding that void galaxies appear bluer, a trend which continues 
down into the most rarefied void regions. Another issue of interest is whether we 
can observe the intricate filigree of substructure in voids, expected as the remaining 
debris of the merging voids and filaments in the hierarchical formation process. 
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Fig. 10. A selection of 10 void galaxies from the SDSS DR7 galaxy redshift survey. These 
galaxies are part of the WSRT HI void galaxy survey. Courtesy: Aragon-Calvo, Platen, 
van de Weygaert. 



5.1. Void galaxy properties 

Void galaxies do appear to have significantly different properties than average field 
galaxies. They appear to reside in a more youthful state of star formation and 
possess larger and less distorted supplies of gas. Analysis of void galaxies in the 
SDSS and 2dFGRS indicate that void galaxies are bluer and have higher specifi c 
star formation rates than galaxies in denser environments M° | ag|o ' | oo | y a | / j | <j | llll | 

While these trends appear to be quite general, controversy persists in the liter- 
ature as to whether or not galaxies in voids differ in their internal properties from 
similar objects in denser regions. Part of this relates to the selection of void galax- 
ies. There are distinct differences between the majority of these galaxies, residing in 
the moderate underdense environment near the boundary, and those in the deepest 
realms, 6 < —0.8, where some 10 % of the void galaxies resides. The most recent 
results on the basis of SDSS DR6 EH anc j 2dFGRS ^ do indicate that these are 
intrinsic differences: at fixed morphology and luminosity the faintest void galaxies 
are intrinsically bluer and have higher star formation rates. It may certainly be true 
that, following the interpretation by ref. because structure grows more slowly in 
the low-density environment of voids the galaxy population in voids is younger. An 
additional effect may be that star formation events may be triggered more easily 
by galaxy encounters in the low velocity dispersion environm ent of v oids ^Sl. 

Interestingly, voids also contain nearly normal ellipticals l ?3 | 19 [ ft seems that 
the red sequence of ellipticals in the color-magnitude diagram from SDSS is almost 
independent of environment con trasting sharply to the significant blueward shift 
of the blue sequence in voids'^. If major mergers produce ellipticals, it is not 
obvious how galaxies in the least dense regions, where interactions are rare, acquire 
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these properties. Careful examination of ellipticals in voids may therefore reveal 
important clues to the balance of nature vs. nurture in galaxy formation. 

5.2. Gas content of void galaxies: the WSRT HI void galaxy 
survey 

Even though a crucial factor for understanding their recent formation and star 
formation history, less is known about the gas content of void galaxies. Possible 
systematic trends in the ratio of hydrogen mass to stellar mass are crucial for 
understanding the influence of the void environment on the galaxies' star formation 
history. The morphology of the gas layer, in particular that of the outskirts, may 
reflect their isolated existence. 

The earlie st sys tematic survey of HI in void galaxies is that of the galaxies in the 
Bootes region HQ2J ft comprised pointed observations of 24 IRAS selected galaxies. 
Sixteen were detected with sufficient amounts of hydrogen, most of which were 
late-type, gas-rich and usually disturbed systems. Interesting was the finding of 18 
new galaxy companions in the vicinity of the known galaxies. Nonetheless, most of 
the observed galaxies were residing in the outer realms of the void and with some 
reason might have been identified with the moderate density environment of walls. 
An important additional HI study of void galaxies involved the search for dwarf 
galaxies in the Local Void ^-A The fact that they were only found at the edge of 




Fig. 11. The SDSS void galaxy selection procedure. Shown is the location of galaxy void078 
in our sample. Top righthand: the SDSS galaxy sky map (z < 0.03), with sky position of 
the galaxy (white dot). The 3-D galaxy sample is translated into a spatial DTFE density 
map, (lower left-hand and centre), and a cosmic spine map (top lefthand). Lower lefthand 
and centre: spatial position of the void galaxy candidate (white dot) with respect to its 
embedding void, seen from two mutually perpendicular directions. The spatial information 
helps to trace the deepest void galaxies. Image courtesy Miguel Aragon-Calvo. 
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Fig. 12. The polar ring galaxy J102819 23+623502.6 (void03). It is one of the most lonely 
galaxies in the Universe, with no NGC and UGC galaxy within 4 Mpc ! Shown are the 
HI intensity map superimposed on SDSS g band optical image (left), the velocity map 
(centre) and the position- velocity diagram along the major axis. The total HI mass in the 
polar ring is 2.0 x 10 9 M Q . See Stanonik et al. 2009a. 



the void and avoid its interior represents a major challenge for our understanding of 
galaxy formation. Some well-known nearby luminous galaxies whose HI layers have 
been studied in detail, ie. NGC6946 and M101, are actually found within the Local 
Void. It is rather intriguing that e.g. NGC6946 is known to undergo a substantial 
level of gas inflow and has a large gas disk bearing the marks of numerous recent 
supernovae El. Galaxies in and around voids are also searched for in the HIPASS 
IS2and ALFALFA survey El. 

The most sensitive and systematic study of gaseous void galaxies is our (ongoing) 
WSRT HI survey of SDSS galaxies. It consists of a sample of around 55 galaxies, 
sufficiently large to probe the morphological diversity of void galaxies and voids of 
different size, shape and underdensity. 

It involves a unique selection procedure in that it is based upon finding void 
galaxies in the deepest underdensities of the DTFE reconstruction of the density 
field ( DTFE: De launay Tessellation Field Estimator) in the (DR7) SDSS redshift 
survey E 4 | 93 | 116 | Fig. |ll| illustrates the void galaxy selection pipeline. Each sample 
galaxy is optically selected from SDSS/DR7 and has to be located in the central 
interior region of a nearby void identified from the SDSS survey. Al sample galaxies 
are at a considerably closer distance than the Bootes void. Most galaxies in our 
sample are blue, a few red (see 10 specimen fig. [T0|). 

Results of the pilot study are reported in refi ll 001 and ref. 11011 With 14 detec tions 
out of of 15 void galaxies, the detection rate of the pilot survey has been 93%HHll 
The one without a detection is the reddest galaxy in the sample. The HI content 
of the 14 detected galaxies ranges from ~ 6.5 x 10 7 M© to ~ 4.8 x 10 9 M Q . Their 
optical radii are small, 5-10 kpc, while they have HI disks that extend far beyond 
the optical and appear to be in regular rotation. Nearly 2/3rd appear asymmetric or 
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slightly warped at the faintest contour levels. Three galaxies turn out to be uniquely 
interesting systems. One is part of a system of three galaxies, all aligned in the plane 
of the sky, stretching out over 57 kpc, and connected by an HI bridge possibly related 
to a tenuous void filament. A second galaxy shows regular disk rotation in its HI 
distribution. However, optically it exhibits a clumpy, knotty morphology similar to 
that seen in chain-galaxies forming at high redshift. It also has a companion. 

The most outstanding specimen of the sample is a polar disk galaxy liffl (see 
fig. IT2|. Amongst the most lonely galaxies in the universe (4 Mpc up to the nearest 
UGC and NGC galaxy), it has a massive, star-poor HI disk that is perpendicular to 
the disk of the central void galaxy. No optical counterpart to the HI disk has been 
found, even though the central galaxy is actively forming stars. The ring reveals a 
regular rotation pattern. 

The galaxy is located within a what looks like a tenuous wall in between two 
large roundish voids. The central optical disk is strongly inclined with respect to the 
surrounding wall (~ 60° — 70°). The unique void environment forces a reappraisal of 
the traditionally assumed formation mechanisms for such polar ring galaxies, by a 
close interaction or merger. Its existence may point to a slow HI accretion dictated 
by the unique environment, possibly at the crossing point of the outflow from the 
two voids. 



5.3. Missing Void Dwarfs: ACDM challenged ? 

With respect to their spatial distribution, void galaxies are found to display a sig- 
nificant luminosity dependence. The consensus is that there is a good statistical 
agreement b etween theory and observation of the large voids defined by L* galaxies 
| 51 | 18 | 46 | 108 | rpjjg observ ed den sity contrast of bright galaxies in large voids, around 
1/10 of the cosmic mean EOlfLU j s consistent with theoretical predictions. 

However, conflicts arise when we look at the distribution and properties of lower 
luminosity/mass objects. From optical and HI surveys, ref.[57l found that the density 
of faint (—18 < Mb < —12) galaxies is only 1/100 that of the mean. This forms a 
marked contrast to the prediction by high-resolution ACDM simulations that voids 
have a densit y of low ma ss (10 9 M o < M < 10 n M o ) halos that is 1/10 that of the 
cosmic mean UHESHnil 

The dearth of dwarf and/or low surface brightness void galaxies represents a 
major challenge for the standard cosmological view of galaxy formation. The sim- 
plest models of galaxy formation predict that voids would be filled with galaxies 
of low luminosity, or galaxies of some other uncommon nature Ref. [78] pointed 
out that the observed salient absence of dwarf galaxies in nearby voids could possi- 
bly involve strong ramifications for the viability of the ACDM cosmology on small 
scales. Branded with the name void-phenomenon, the issue has been marked as a 
"a crisis for the CDM model" since Peebles strongly emphasized the notion that 
nearby voids do NO T appear to contain the expected population of dwarf and/or 
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low surface brightness galaxies. 

Various processes and influences have been forwarded as an explanation for the 
missing void dwarfs within the context of existing galaxy formation theories. These 
involve environmental properties of dark matter halos and astrophysical procesess, 
ranging from supplementary baryonic physical processes and radiation feedback 
processes. 

The nature of the environmental influence on the formation of halos has been 
addressed in a variety of studies. High- resolution simulations now clearly predict 
dependence of the s hape, concentratio n, spin, and mean age of formation of halos 
on their environment BHEHHH15II11B particularly so for lower mass halos. Because 
the mass function of halos shifts to low er masses in voids these effects should be most 
noticeable among void galaxies | 34 | 33 |_ Ref . [301 and ref . 1 1 081 modelled the dependence 
of the properties of galaxies on the halos they inhibit to infer the properties of voids 
in the galaxy distribution. According to ref. I108[ halo models would explain away 
the void phenomenon. However, the implied spatial segregation between galaxies of 
different luminosities in and near voids does not bear much resemblance to reality. 
The issue still seems far from solved and progress will depend largely on insights 
offered by new observational data. 

Amongst the most plausible astrophysical explanations for the absence of dwarfs 
in voids is the ionizing influence of the UV background in suppressing the star 
forming capacity of low mass halos HI However, ref. US] showed that these effects 
would not explain the distribution of dwarf galaxies in voids. Also the suggestion 
that supernova feedback would be instrumental in evicting gas from the shallow 
potential wells of the dwarfs does not seem to offer an explanation. Simulations by 
ref. 11031 showed that most internal properties of dwarf galaxies can be reproduced 
without invoking supernova feedback. 
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